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Abstract. Static spherically symmetric anisotropic source has been studied for the Einstein- 
Maxwell field equations assuming the erstwhile cosmological constant A to be a space-variable 
scalar, viz., A = A(r). Two cases are examined out of which one reduces to isotropic sphere. 
The solutions thus obtained are shown to be electromagnetic in origin as a particular case. It is 
also shown that the generally used pure charge condition, viz., p + p r = is not always required 
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£T) ' for constructing electromagnetic mass models. 
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The cosmological constant A, related to the energy of space, introduced by Einstein in general 
relativity has become of late very significant from the view point of cosmology. Though Einstein 
ultimately abandoned it stating that it was a "blunder" in his life but Tolman keeps it as a con- 
stant quantity in his field equations even in 1939's famous work related to astrophysical system 
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1. Introduction 



(Tolman, 1939). It is also, in favor of keeping A, argued by Peebles and Ratra (2002) that like all 
energy, the zero-point energy related to space has to contribute to the source term in Einstein's 
gravitational field equations. However, it is being gradually felt that the erstwhile cosmological 
constant A is indeed a scalar variable dependent on time rather than a constant as was being 
believed earlier. Recently, this variation in cosmological constant is also observationally con- 
firmed due to the evidence of high redshift Type la supernova (Perlmutter et al, 1998; Riess et 
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al., 1998) for a small decreasing value which is < 10~ 56 cm~ 2 at the present epoch. Obviously, 
once the A becomes a scalar its dependence need not be limited only to time coordinate (as in 
cosmology). Since it enters in the field equations as a variable, it must be dependent on space 
coordinates as well. Therefore, in general, the A is a scalar variable dependent on, either or 
both, space and time coordinates. We believe that just as in cosmology the dependence of A on 
time has been found to be of vital importance playing a significant role now, its dependence on 
space coordinates is equally important for astrophysical problems, in particular, the problems 
of small dimensions like that of an extended electron. With this view point, we consider here a 
charged anisotropic static spherically symmetric fluid source of physical radius, a, by introduc- 
ing a scalar variable A dependent on the radial coordinate. The field equations thus obtained, 
under certain mathematical assumptions, yield a set of solutions which has a kind of another 
historical importance, known in the name of Electromagnetic Mass Models (EMMM) in the lit- 
erature (Feynman et al, 1964). The effective gravitational mass of these models depends on the 
electromagnetic field alone (the effective gravitational mass vanishes when the charge density 
vanishes). Such models have been studied by several authors (Tiwari et al. 1984, 1986, 1991; 
Gautreau, 1985; Gr0n, 1985, 1986a, 1986b; Ponce de Leon, 1987a, 1987b, 1988; Tiwari and Ray, 
1991a, 1991b, 1997; Ray et al, 1993; Ray and Ray, 1993). All these EMMM, however, have 
been obtained under a special assumption p + p = 0, where p is the matter-energy density and 
p is the fluid pressure under the general condition that p > and p < 0. This type of equation 
of state implies that the matter distribution is in tension and hence the matter is known, in 
the literature, as a 'false vacuum' or 'degenerate vacuum' or 'p-vacuum' (Davies, 1984; Blome 
and Priester, 1984; Hogan, 1984; Kaiser and Stebbins, 1984). A natural question arises whether 
there exists any EMMM where this condition is violated, i.e., when p + p / 0. This is the main 
motivation of the present investigation and here we have shown that even for p + p / EMMM 
can be constructed. However, to carry out this plan we have assumed here that the so-called 
cosmological constant is dependent on the space coordinate (Chen and Wu, 1990; Narlikar et 
al, 1991; Ray and Ray, 1993; Tiwari and Ray, 1996). This will help us to have several class of 
solutions related to EMMM. 

2. Einstein-Maxwell field equations 

Let us consider a spherically symmetric line element 

ds 2 = g ij dx i dx j = e v{r) dt 2 - e x{r) dr 2 - r 2 (d6 2 + sin 2 8d<f> 2 ), = 0, 1,2,3). (1) 
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Now, the Einstein field equations for the case of charged anisotropic source are 

G*i = W 3 - \g\R + g^A = -«pT*/ m) + r, (em) ], (2) 

where T 1 ^ and T l< f m ^ are respectively the energy-momentum tensor components for the 
anisotropic matter source and the electromagnetic field and are given by 

r/ m) = (p + Pl _)u l Uj - VL tf 3 + (p ± - Pr )rf m (3) 

with UiU 1 = —r]iif = 1, and 

Tl (em) = _L { _ FjkF ik + ±-tf F kl F kl ). (4) 

The Maxwell electromagnetic field equations are given by 

[(-g^F^^An-ri-gf 2 , (5) 

and 

F[i jM = 0, (6) 

where the electromagnetic field tensor Fij is related to the electromagnetic potentials through 
Fij = Aij — Aj t i which, obviously, is equivalent to the equation (5), viz., F^j^] = 0. Further, 
u % is the 4- velocity of a fluid element, J 1 is the 4-current satisfying J 1 = au % , where a is the 
charge density, and k = 8n (in relativistic unit G = C = 1). Here and in what follows a comma 
denotes the partial derivative with respect to the coordinates (involving the index). 
The Einstein-Maxwell field equations (2)-(6) corresponding to anisotropic charged source with 
cosmological variable, are then given by 

e- x (X'/r - 1/r 2 ) + 1/r 2 = 87rT° = 87rp + E 2 + A, (7) 

er x (u'/r + 1/r 2 ) - 1/r 2 = -8t:T 1 1 = 8vrp r - E 2 - A, (8) 

e" V + ^' 2 /4 - ^A'/4 + [y' - A')/2r] = -8ttT 2 2 = -8irT 3 3 = 8ir P± + E 2 - A, (9) 

{r 2 E)' = 4nr 2 ae x / 2 . (10) 
The equation (10) can equivalently, in terms of the electric charge q, be expressed as 

q(r)=r 2 E{r)= [ Airr 2 ae x / 2 dr (11) 
Jo 
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where p r , p± and E are the matter-energy density, radial and tangential pressures and intensity 
of the electric field respectively. Here prime denotes derivative with respect to radial coordinate 
r only. 

The equation of continuity, T l j; i = for the field equations (2) - (6), is given by 

4- \p r - (E 2 + A)/8vrl + (p + Pr y/2 = E 2 jlvr + 2(p ± - p r )/r. (12) 
dr l J 

Now, we assume the relation between the radial and tangential pressure as 

p ± = np r , (n ^ 1). (13) 

Assuming further that the radial stress T 1 1 = (Florides, 1987; Kofinti, 1985; Gr0n, 1986b; 
Ponce de Leon, 1987b; Tiwari and Ray, 1996), one gets 

y' = (e A - i)/r, (14) 

p r = (E 2 + A)/8?r. (15) 
Using equations (13) - (15), in equation (12), we get 

p+p r = [(n+l) J E 2 + (n-l)A)]/27r(e A -l). (16) 

Similarly, equations (7) and (8) yield, 

e- x {u' + \') = 8irr(p + Pr ). (17) 

Again, equation (7) together with equation (15), gives 

e ~ x = l-2m(r)/r, (18) 

where m(r), called the gravitational mass, takes the form 

rr 

m(r) = M(r) + p(r) = Air [p + p r ]r 2 dr, (19) 

Jo 

the Schwarzschild mass and the mass equivalence of electromagnetic field, respectively, being 
defined as 

rr rr 

M(r) = 47r / pr 2 dr, p(r) = An p r r 2 dr. (20) 

JO JO 

Now, from the above equation (19) it is easily observed that the condition p + p r = 0, yields a 
flat space-time through the equations (18) and (14) and has been considered by us in another 
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context (Tiwari et al., 2000). Hence, the non-trivial solutions exits here for the case p + p r / 
only. 

3. Solutions for the static charged fluid spheres 

Let us now solve the equation (16) under the constraint p + p r ^ 0, assuming different mathe- 
matical conditions. As the equation (16) is involved with two physical parameters so unless we 
specify one parameter it is not possible to solve the equation (19). Let us, therefore, assume the 
following two cases which will yield solutions with physically interesting features as the analysis 
in Section 4 demonstrates. 
Case I: 

We now make the choice 

A = E 2 - NA , (21) 

where N is an integer and Ao is the erstwhile non-zero cosmological constant. 
With the help of equations (16) and (21), the equation (19) takes the form 

m(r) = 2 [[2nE 2 - (n - l)NA ]r 2 dr/(e x - 1). (22) 
Jo 

To make equation (22) integrable we assume that 

E 2 = q 2 /^ = [k(e x - 1)(1 - R 2 ) + (n - l)iVA ]/2n, (23) 

where A; is a constant and R = r/a, a being the radius of the sphere. This particular choice for 
the electric intensity generates a model for charged sphere which is physically very interesting 
as it is related to EMMM as will be seen later on. Thus, the solution set is given by 

e~ x = 1 - AR 2 (5-3R 2 ), (24) 

e u = (1 - 2A) 5 / 4 e A/4 exp[5 J Btan- 1 S(6 J R 2 - 5) - tan- l B/2], (25) 

p r = p± / n = [k(e x - 1)(1 - R 2 ) - A^Aol/Svrn, (26) 

p = [k(l + 4n - e A )(l - R 2 ) + iYA ]/8vrn, (27) 

where 

A = 4ka 2 /15, B = [A/ (12 -25 A)} 1 / 2 . (28) 
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Now, the exterior field of a spherically symmetric static charged fluid distribution described by 
the metric (1) is the unique Reissner-Nordstrom solution given by 



2m q 2 




2m q 2 




dt 2 - 











dr 2 -r 2 {d6 2 + sin 2 ed<t> 2 



(29) 



Then, application of the matching condition on the boundary r = a yields the total effective 
gravitational mass which is given by 



, , 4, , q(a) 2 



(30) 



Hence, in terms of the total gravitational mass, the total electric charge and the radius of the 
sphere the constants k, A and B can be expressed as 



and 



k = l5(2am-q 2 )/8a i , 
A = (2am - q 2 )/2a 2 , 

B 2 = (2am - q 2 )/[2Aa 2 - 25(2am - q 2 



(31) 
(32) 

(33) 



Comparing the equation (19) (and hence via the equation (20)) with the above equation (30) we 
can easily recognize the first term (4/ca 3 /15) in the right hand side as the total Schwarzschild 
mass whereas the second term (q 2 /2a) is the total mass equivalence of the electromagnetic field. 
Now, k, in the first term of equation (30), is implicitly related to the charge q as is evident from 
the equation (23). Therefore, in principle, the gravitational mass m is of purely electromagnetic 
origin. 
Case II: 
Here the choice is 

A = iVA - E 2 . (34) 
Then, from the equations (16) and (34), the gravitational mass (equation (19)) reduces to 

m(r) = 2 ( r [2E 2 + (n - l)A^A ]r 2 dr/(e A - 1). (35) 
Jo 

For the further assumption 

E 2 = q 2 /r 4 = [k(e x - 1)(1 - R 2 ) - (n - l)iVA ]/2, (36) 
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we have the solution set 

e~ x = 1-AR 2 (5-3R 2 ), (37) 

e u = (1 - 2A) 5/4 e x/4 exp[5Btan' 1 B(6R 2 - 5) - tan~ 1 B/2], (38) 

p = NA /8tt, (39) 

p=[4k(l-R 2 )-NA ]/87T. (40) 

We see that A and v retain the same form as in case I and hence the total gravitational mass 
is also given by the equation (30). Further, it is observed that the present case automatically 
reduces to an isotropic one as the pressure p does not associated with the anisotropic factor n. 

4. Physical properties of the static charged fluid spheres 
Case I: 

(1) The equation (23) indicates that for getting a direct dependence of k upon q one can admit 
the following relaxation such that (i) N = when Aq ^ and (n — 1) / 0, (ii) (n — 1) = when 
N ^ and Aq ^ and (iii) Ao = when (n — 1) / and N ^ 0. The third possibility seems 
to contradict the observational results related to the Supernova type la where the cosmological 
constant is found to be a non-zero value whereas the second one provides an isotropic case. 
Then, suitably opting for the sub-case (i), viz., N = one obtains 

q 2 = kr 4 (e x -l)(l-R 2 )/2n, (41) 

Pr =p±/n = k{e x - 1)(1 - i? 2 )/8vrn, (42) 

p = k(l+4n-e x )(l-R 2 )/8irn. (43) 

Thus, for vanishing electric charge the gravitational mass in the equation (30), including all the 
physical parameters (viz., pressures and density), vanishes and one obtains EMMM. 

(2) The central and the boundary pressures are found to be equal here i.e., p r (0) = p±(0)/n = 
—NAo/8im and p r (0) = p±(0)/n = —NAo/8im respectively whereas the respective densities are 
{Ank + N Aq) / 8im and NAq^ttu. So, the present model has a constant pressure throughout the 
sphere though the density decreases from center to boundary. For the value N = 0, however, 
we have zero pressure, both at the center and boundary, and the density decreases from the 
central value k/2iT to zero at the boundary. Thus, with N = the present model goes to a 
physically well-behaved static charged dust case. However, we would like to investigate how the 
pressure and the density behave in between the center to boundary. To see this, as an example, 
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we have plotted the graph for the extended classical electron of Lorentz type (1904) in Figure 
^ It shows the nature of variation of the fluid pressure and the energy density with radius from 
the center of the matter distribution to the boundary. We observe that though the pressure 
is zero at the center as well as at the boundary but the in-between feature is very interesting. 
Starting from zero it increasing slowly and then after attaining a maximum value drops steadily 
to zero. In the case of density the curve starts from the non-zero central value k/2~K and then 
smoothly decreases to zero at the boundary. However, it is to be noted here that for the electron 
with radius a = 10 -16 cm, mass m = 6.76 x 10~ 56 cm and charge q = 1.38 x 10 -34 the value 
of the constant k becomes a negative quantity with —3.568215 x 10~ 4 . This particular aspect 
one should keep in mind for the explanation of the internal structure of classical electron with 
respect to the fluid pressure and the energy density. Thus we see that the density starts from the 
down of the axis due to this negative value attached to it. In the Figure |2] the pressure-density 
profile is shown for different values of the anisotropic factor n. 




(3) In the above analysis for the general value of N we have seen that the fluid pressure is 
altogether negative whereas the density is a positive quantity. Now, it can be observed from the 
equations (26) and (27) that p+p r / except at the boundary where it is equal to p r (a) = —p(a), 
where p(a) > and p{a) < as seen earlier. Otherwise, it will have a general value fc(l — i? 2 )/27r 
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Figure 2: The fluid pressure and the energy density profile is shown for a classical electron with 
radius a = 10~ 16 cm, mass m = 6.76 x 1CP 56 cm and charge q = 1.38 x 10~ 34 cm. 

. The central value is then p r (0) = —p(0) where p(0) = p(0) — k/2ir. As the model demands for 
p > and p < so the condition to be satisfied here is p(0) > k/2ir. The general condition for 
the negative pressure and positive density is then p > k(l — R 2 )/2tt for all r < a. These results 
are also true for the sub-case N = 0. 



Case II: 

(1) Here for N = the solution set regarding the pressure and density (vide equations (39) and 
(40)) reduces to 

V = 0, (44) 
p = fc(l - i? 2 )/27r, (45) 

when the electric charge is given by 

q 2 = k(e X - 1)(1 -i? 2 )r 4 /2. (46) 
Thus, as in the previous case, for q = we get k = which in turn makes mass and density 
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everything zero and also the space-time becomes flat. Thus, the model presented here is an 
EMMM. 

(2) In the present case also the central and the boundary pressures are equal with a value 
NAq/8tt and the respective densities are (4A; — NAq)/8tt and — iVAo/87r. The pressures become 
zero for N = 0, both at the center and boundary, and densities have the positive central value 
k/2ir whereas the boundary-value is zero. Thus, we again get a physically interesting charged 
dust case with N = 0, only difference with the previous Case I is that this is now a case of 
isotropic fluid sphere. Figure |3] shows the nature of variation of energy density with radius from 
the center of the matter distribution to the boundary for an extended classical electron. Here 
also the behavior is regular and well defined. 




0.2 0.4 0.6 0.8 1 

radius (cm) xio" 16 



Figure 3: The energy density as a function of radius is plotted for a classical electron with radius 
a = 10 -16 cm, mass m = 6.76 x 10 -56 cm and charge q = 1.38 x 10~ 34 . 

(3) In the present case also, by virtue of equations (39) and (40), p + p T ^ which reads here 
as p + p r = k{\ — R 2 )/2tt. Here the central value, at r = 0, is p = —p where p = (p — k/2ir) 
and the boundary one is p = —p. Due to negativity of the density here the condition on the 
pressure to be imposed is p > k/2n. Thus, the present Case II clearly provides an EMMM even 
with a positive pressure and therefore contradicts the comment made by Ivanov (2002) that 
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"... electromagnetic mass models all seem to have negative pressure." The same result, i.e. the 
positivity of pressures are also available in some cases of the work done by Ray and Das (2002) 
related to EMMM. However, the explanation given here is valid for any positive value of N and 
so the situation may completely be opposite with any negative value of N. At this stage, we feel, 
we should not put any restriction on the choice of the value of N. This is because, in general, for 
a fluid sphere we should have p > and p > so that the weak energy conditions are satisfied. 
But there are also some special situations available within the spherical system (particularly, 
in the case of electron with the radius ~ 10~ 16 cm) where the energy condition is violated 
due to negative energy density (Cooperstock and Rosen, 1989; Bonnor and Cooperstock, 1989). 
Thus, choosing the proper signatures of ./V we can have a class of models with diverse characters. 

5. Role of A: Previous and Present Status 

The cosmological constant was introduced by Einstein in his field equation to obtain a static 
cosmological solution because of the fact that due to gravitational pull everything will collapse 
to a point and hence a un-wanting situation of singularity will take place. However, he was not 
satisfied with this new physical quantity as it seemed to violate Machian principle which he tried 
to incorporate in the framework of his general theory of relativity. He thus, ultimately rejected 
it mainly for two reasons: (i) that the theoretical work of de Sitter showing that the Einstein's 
field equations admitted a solution for empty Universe and (ii) that the experimental discovery 
of expanding Universe by Hubble. 

As stated in the introduction, the concept of cosmological constant has been revived recently in 
the case of early Universe scenario and even in particle physics. It is gradually being felt that 
A, the erstwhile cosmological constant is available rather than a constant, as was being believed 
earlier, varying with space or time or both (Ray and Ray, 1993; Tiwari and Ray, 1996; Tiwari 
ct al., 2000). Further, A may be positive or negative (by imposing the condition that its value is 
not equal to zero). For instance, according to Zel'dovich (1968) the effective gravitational mass 
density of the polarized vacuum is negative. Similarly, the equation of state p+p = 0, employed 
by Tiwari et al. (1984) to construct EMMM as a solution of Einstein-Maxwell field equations, 
provides negative pressure. It may be emphasized here that positive density has significant, 
rather major role in inflationary cosmology whereas negative density has influence on elemen- 
tary particle models. The gravitational mass inside the spherical charged body is negative for 
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r < 5a/4, where r is the radial coordinate and a is the radius of the sphere (Gautreau, 1985). It 
is argued by Gr0n (1986a, b) that this negative mass and the associated gravitational repulsion is 
due to the strain of the vacuum because of vacuum polarization. He also argue that if a vacuum 
has a vanishing energy, then its gravitational mass will be negative and the observed expansion 
of the universe may be explained as a result of repulsive gravitation. Now, if we consider a 
negative A having a repulsive nature as was considered by Einstein then this gets the same 
status of negative pressure and also can be identified with the Poincare' stress. This repulsive 
gravitation associated with negative A can also be explained as the source of gravitational blue 
shift (Gr0n, 1986a). On the contrary, positive A will be related to gravitational red shift. It 
may be also pointed out that according to Ipser and Sikivie (1984) domain walls are sources of 
repulsive gravitation and a spherical domain wall will collapse. To overcome this situation the 
charged "bubbles" with negative mass keep the wall static and hence in equilibrium. In this 
regard, we may also add that A, via repulsive gravitation, is related to domain walls and playing 
an important physical role. 

Very recent observations conducted by the SCP and HZT (Perlmutter et al, 1998; Riess et al, 
1998; Filippenko, 2001; Kastor and Traschen, 2002) show that the present value of A is positive 
one and hence related to the repulsive pressure. It is believed that the present state of accelera- 
tion dominated universe is due to the driven force of this A. It is to be noted that the negative 
A corresponds to a collapsing situation of the universe (Cardenas et al, 2002). 

6. Conclusions 

(i) In both the above cases I and II, it is possible to show that EMMM also can be obtained, in 
principle, using the constraint p + p r ^ 0. This particular point remained unnoticed by Gr0n 
(1986a, b) and Ponce de Leon (1987a, b), both. However, as p + p r = is related to vacuum 
polarization, black hole physics and inflationary cosmology (Gr0n, 1985, 1986a; Wenda and Shi- 
tong, 1985a, b; Guth, 1981; Linde, 1984), so the present model , in general, is in contradiction 
to those phenomenological explanations. 

(ii) It can be noted that in terms of energy-momentum tensor of the fluid the condition p+p r = 
implies T 1 ! = T°o (Herrera and Varela, 1996) whereas p + p r ^ constraint may be expressed 
as T 1 ! = as we have adopted in the present approach. It is also interesting to note that 
p + p r = and hence T l \ = T° can be expressed in terms of the metric tensors (vide equation 
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(1)) as gooffii = ~~ 1- A coordinate-independent statement of this relation is obtained by Tiwari 
et al. (1984) by using the eigen values of the Einstein tensor G l j. 

(hi) Another point is that in both the cases we have taken the assumptions in a way so that 
the cosmological variable A does not vanish rather may be at most equal to Ao, the erstwhile 
cosmological constant having a finite non-zero value. Without considering Ao we will have A = 
at the boundary r = a which is a bit unphysical and may create difficulties, such as to entropy 
like problems (Beesham, 1993). In this regard, one can observe that the solutions obtained by 
Gr0n(1986a, b) and Ponce de Leon (1987a, b) represent a neutral system, viz., though the net 
charge is not zero but the charge on the surface of the spherical system vanishes. The models 
of the present paper, in general, do not correspond to this situation because of the fact that 
the cosmological parameter A does not vanish at the boundary, rather, its value is Ao at r = a. 
Therefore, the present solutions correspond to charged sphere. Of course, for N = 0, like Gr0n 
(1986, b) and Ponce de Leon (1987a, b), we have neutral spheres (vide equation (41) of case I 
and equation (46) of case II). Thus, we have a class of solutions related to charged as well as 
neutral systems depending on the values of N. 

(iv) We have, as a special case, applied the present models to the classical electron of Lorentz 
type. However, it seems that a more realistic application of the models with variable cosmologi- 
cal constant is possible in the case of charged massive astrophysical systems, like neutron stars. 
This aspect will be carried out in detail in the future investigations. 
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